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Introduction 

1,2-Dibromocthane  or  ethylene  dibromide  (EDB)  is  a  common 
groundwater  and  soil  pollutant.  EDB  was  formally  used  as  an 
antiknock  additive  in  leaded  gasoline  and  aviation  fuel.  It  was 
also  used  as  a  pesticide  for  citrus  and  grains  until  its  use  as  a 
fumigant  was  banned  by  the  U.S.  Environmental  Protection 
Agency  (EPA)  in  1984.  Contaminated  groundwater  sites  are  usu¬ 
ally  associated  with  pesticide  use  or  fuel  spills.  U.S.  EPA  has 
established  a  maximum  contaminant  level  (MCL)  of  50  parts  per 
trillion  (ppt)  for  EDB  in  drinking  water,  due  to  its  toxicity  and 
carcinogenicity.  The  State  of  Massachusetts  has  set  a  MCL  of  20 

ppt. 

Zero  valent  zinc  was  first  reported  as  a  reducing  agent  in  the 
dihalo-elimination  of  EDB  in  1874  (Baciocchi  1983).  More  re¬ 
cently,  Rajagopal  and  Burris  (1999)  found  that  when  EDB  in 
water  was  exposed  to  Fe^  the  only  observed  products  were  ethyl¬ 
ene  and  bromide 

BrH2CCH2Br+Fc^^^H2CCH2  +  2Br+Fe2  (1) 

Recently  zero  valent  iron  (ZVI)  has  been  applied  to  the  reme¬ 
diation  of  contaminated  groundwater  (Gillhiam  and  O’Hannesin 
1994).  A  great  deal  of  work  has  been  done  on  the  reduction  of 
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chlorinated  organics,  nitrogenous  organics,  and  dissolved  metals 
by  ZVI  (Matheson  and  Tratnyek  1994;  Agrawal  and  Tratnyek 
1996;  Siantar  et  al.  1996;  Hundal  et  al.  1997;  Sayles  et  al.  1997; 
Shokes  and  Moller  1999).  Permeable  ZVI  reactive  barriers  have 
been  put  in  place  in  a  number  of  locations  (O’Hannesin  and  Gill- 
ham  1998;  Puls  et  al.  1998;  Gavasker  et  al.  1999).  A  permeable 
reactive  barrier  (PRB)  is  a  subsurface  media  (catalyst,  reactive 
surface,  or  biological  system)  of  high  porosity  that  intercepts  a 
plume  of  contaminated  groundwater  in  situ  and  removes  the  con¬ 
taminant  of  interest.  The  reduction  rates  of  chlorinated  ethencs 
[e.g.,  trichloroethylene  (TCE)  and  tetrachloroethylene  (PCE)] 
with  ZVI  are  relatively  slow,  however,  ZVI  based  PRB’s  can  be 
used  because  low  groundwater  flow  rates  provide  sufficiently 
long  residence  times  for  complete  dechlorination. 

EDB  reacts  faster  than  many  chlorinated  ethylencs.  In  similar 
experiments  (5g  Fisher  40  mesh  iron  and  100  mLDI)  the  half  life 
of  EDB  was  found  by  Rajagopal  and  Burris  (1999)  to  be  L3  h, 
and  the  approximate  half  lives  of  1,1,1-trichloroethane,  TCE  and 
PCE  were  1.6,  24,  and  120  h,  re.spcctively  (Johnson  et  al.  1996; 
Campbell  et  al.  1997).  The  faster  reaction  kinetics  of  EDB  reduc¬ 
tion  by  ZVI  would  allow  for  shorter  residence  times,  possibly 
expanding  the  engineered  treatment  options  available  for  EDB 
contaminated  groundwaters. 

Most  in  situ  applications  of  ZVI  PRB’s  require  excavation  of 
aquifer  material  and  replaccntent  with  iron  filings,  which  have 
cost  and  logistical  limitations.  If  a  relatively  short  reaction  zone  is 
adequate  for  groundwater  flow  rates,  the  use  of  hydraulic  fractur¬ 
ing  to  emplace  the  ZVI  PRB  (Hocking  et  al.  1998)  may  be  fea¬ 
sible. 

Above  ground  packed-bed  ZVI  reactors  may  accommodate  the 
higher  flow  rates  required  for  the  more  traditional  pump-and-treat 
operations.  In  the  above  ground  treatment  system  the  iron  could 
be  more  easily  replaced  as  needed  due  to  loss  of  reactivity  or 
clogging.  Additionally,  the  syslcm  could  be  relocated  as  neces¬ 
sary. 

In  this  study,  the  kinetics  and  mass  transfer  characteristics  of 
EDB  reduction  in  laboratory  and  pilot  scale  ZVI  packed-bed  re¬ 
actors  were  examined.  A  model  incorporating  reaction  kinetics 
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and  mass  transfer  limitations  was  developed  and  used  to  analyze 
the  results  from  the  experimental  systems.  The  feasibility  of  treat¬ 
ing  EDB  contaminated  groundwater  using  (1)  an  above  ground 
treatment  system,  and  (2)  a  naiTOw  width  in  situ  reactive  wall 
(that  could  be  emplaced  by  hydraulic  fracturing  or  jetting)  were 
evaluated  using  the  model. 

Materials  and  Methods 

Chemicals 

The  cast  iron  (GX-27,  Master  Builders,  Cleveland)  was  used  as 
received.  The  cast  iron  filings  were  nonuniform  flakes  and  cylin¬ 
ders  with  an  average  particle  diameter  dp  of  0.36  mm.  The  sur¬ 
face  area  (SA)  measured  using  the  nitrogen  BET  method  (Micro- 
metrics,  Flowsorb  2300,  Norcross,  Ga.)  was  0.79  mVg.  The  sand 
(Ottawa  Sand)  was  water  rinsed  and  oven  dried. 

Water  used  in  the  packed-bed  reactor  studies  was  groundwater 
contaminated  with  EDB  collected  from  the  FS-12  treatment  site  at 
the  Massachusetts  Military  Reservation,  Mass.  (FS-12  MMR). 
The  pH  ranged  from  6.3  to  7.3  and  dissolved  oxygen  concentra¬ 
tions  ranged  from  6.3  to  8.7  mg/L.  The  influent  water  temperature 
varied  seasonally  from  10.2  to  12.0^C.  On-going  monitoring  of 
the  FS-12  site  gave  the  following  water  quality  parameters; 
total  dissolved  solids=  46.4  mg/L,  total  suspended  solids=2.3 
mg/L,  total  organic  carbon =0.3 3  mg/L,  Fe= 40.81  p-g/L,  Mn=16 
|jL,g/L,  oxidation-reduction  potential =260  mV. 

Dibromoethane  (Photrex  grade)  was  obtained  from  J.  T.  Baker. 
Acetonitrile  [high-pressure  liquid  chromatography  (HPLC) 
grade],  Methanol  (HPLC  grade),  and  p-dichlorobenzene  (DCB, 
reagent  grade)  were  obtained  from  Fisher  Scientific.  Pentane  [tri- 
halomethane  (THM)  grade]  was  obtained  from  Aldrich. 

Batch 

The  batch  kinetic  experiments  were  conducted  in  160  mL  crimp 
top  glass  vials  sealed  with  Teflon  faced  rubber  septa,  100  mL  of 
unbuffered  de-ionized  (DI)  water,  anaerobic  headspace,  and  5  g 
iron  were  spiked  with  1-10  |jlLEDB  (0.011-0.12  mM).  The  vials 
were  rotated  at  24  rpm  (0.51  m/s)  in  the  dark.  Reaction  tempera¬ 
tures  were  maintained  in  incubators,  at  5,  10,  20,  30,  and  40 
(±1)°C.  The  mass  balance  experiments  were  performed  at  20X. 
Headspace  samples  (100  ixL)  were  periodically  obtained  using  a 
gas-tight  syringe  and  analyzed  by  GC/FID.  Control  samples  with 
no  iron  were  run  concurrently  and  showed  negligible  loss  of  EDB 
over  several  days.  All  experiments  were  run  in  duplicate. 

Sorption  experiments  were  run  in  zero  headspace  batch  sys¬ 
tems  with  duplicate,  sacrificial  samples.  Iron  filings  (1.5  g)  were 
placed  in  15  mL  serum  vials,  vials  were  filled  completely  with  DI 
water  (^14.8  mL)  in  an  anaerobic  chamber  (H2  5%,  N2  95%)  and 
crimp  sealed  with  Teflon  faced  rubber  septa.  Immediately  after 
being  sealed  the  vials  were  removed  from  the  chamber  and  spiked 
with  methanol  solutions  of  dibromoethane  ('^1,260  mg/mL)  to 
give  a  final  concentration  of  ~85  mg/L  (0.45  mM).  The  water 
volume  of  each  vial  was  determined  gravimetrically.  Control 
samples  contained  no  iron.  The  samples  were  mixed  by  axial 
rotation  on  a  roller  drum  at  24  rpm,  at  20±1°C  in  the  dark. 
Duplicate  control  and  experimental  samples  were  removed  after 
varying  periods  of  time  and  analyzed  to  determine  aqueous  and 
total  EDB  concentrations.  Duplicate  20  |jlL  aliquots  of  the  aque¬ 
ous  phase  were  sampled  from  each  vial  and  spiked  into  1.0  mL  of 
acetonitrile  with  10  mg/L  /7-di chlorobenzene  as  an  internal  stan¬ 
dard  and  vortexed.  Total  system  EDB  concentrations  were  deter¬ 


mined  by  transferring  the  aqueous  content  of  each  vial  into  a  40 
mL  vial  containing  10  mL  of  acetonitrile  via  cannule.  The  remain¬ 
ing  solids  were  washed  with  two  2.5  mL  portions  of  acetonitrile 
(spiked  with  internal  standard)  and  the  washing  solvents  added  to 
the  40  mL  vial.  A  100  jxL  aliquot  was  spiked  into  0.90  mL  of 
acetonitrile  and  analyzed  via  GC/ECD. 

The  difference  between  the  total  and  aqueous  phase  EDB 
masses  was  assumed  to  be  the  EDB  mass  bound  to  the  metal 
surface.  The  sorption  isotherm  (20°C)  was  made  by  plotting  the 
bound  concentration  (|xg/m^)  versus  the  aqueous  concentration 
(mol/m^)  as  the  reaction  progressed. 

Laboratory  Scale  Packed-Bed  Reactor 

The  lab  scale  reactor  (Fig.  1)  was  a  column  made  of  Plexiglas 
[19.7  cm  lengthX  2.54  cm  inner  diameter  (i.d.)]  with  stainless 
steel  sampling  ports  and  fittings.  Stainless  steel  screens  were 
placed  at  the  column  inlet  and  outlet  followed  by  22  g  sand,  254 
g  iron,  and  28  g  sand  (inlet  end  down).  Sampling  ports  1  and  10 
were  located  in  the  sand  zones.  Sampling  ports  2-9  were  located 
in  the  iron  zone  (12.7  cm  in  length).  The  sampling  ports  were 
0.158  outer  diameterX  0.127  cm  i.d.  tubing  with  vertical  half 
slots,  attached  to  the  column  using  0. 1 58  x  0. 1 58  cm  bore  through 
male  NPT  Swagelok  fittings  with  Teflon  inserts  to  minimize  dead 
volume.  Mininert  valves  were  placed  on  the  end  of  sampling  tube 
to  facilitate  sampling. 

The  pore  volume  was  determined  by  pumping  a  known  vol¬ 
ume  of  water  through  the  dry  columns  and  collecting  the  effluent, 
the  difference  was  taken  to  be  the  total  column  pore  volume 
(  I86±3  mL).  As  H2  was  produced  during  operation  of  the  iron 
column,  pore  water  was  displaced;  tiius  the  actual  pore  volume 
during  operation  was  slightly  less  than  ihai  determined.  The  po- 


86 /JOURNAL  OF  ENVIRONMENTAL  ENGINEERING /JANUARY  2002 


cv 


Fig.  2.  Schematic  of  zero  valent  iron  pilot  scale  packed-bed  reactor.  HV=hand  valve,  SV^^solenoid  valve,  FI=^Rotoineter,  Pump-variable  speed 
pump,  Airtrap^ water  hammer  suppressor,  PRV=pressurc  relief  valve,  PT= pressure  gauge,  PT^electronic  pressure  sensor,  INF  SP=influcnt 
sample  port.  Flow  sensor,  Sacrihcial  DO  Scrubbers^one  in-line,  one  in  reserve,  packed-bed  canisters  (11),  SP=Sample  ports  (12), 
DP^diffcrcntial  pressure  gauge,  bypass  valve,  and  CV -check  valve.  _ _ _ _ 


rosity  (e)  was  found  by  dividing  the  pore  volume  by  the  empty 
bed  volume  of  the  column  (e  =  0.53 ±0.05).  This  number  in¬ 
cludes  the  porosity  in  the  sand  layers,  because  of  the  relatively 
small  amount  of  sand  used,  the  error  in  the  mea.sured  porosity  was 
assumed  to  be  small. 

The  column  was  sampled  at  various  time  intervals  at  sampling 
ports  I -10.  Samples  were  drawn  from  the  top  sampling  port  first, 
sequentially  to  the  influent.  Samples  were  collected  using  a  150 
|ulL  syringe  and  the  first  100  |jlL  was  discarded  (dead  volume  ^^25 
|xL).  The  50  p,L  sample  was  then  diluted  as  described  below  and 
analyzed  by  GC/ECD. 

Flow  was  maintained  with  either  an  ISCO  2350  or  Waters  510 
(Waters,  Milford  Mass.)  HPLC  pump,  at  rates  of  1.00±0.05, 
5.00±0.25,  and  10.0±0.5  mL/min;  flow  rates  w'ere  measured  pe¬ 
riodically  through  out  the  experiment  to  ensure  consistency.  The 
influent  water  was  groundwater  collected  at  FS-12  MMR,  spiked 
with  stock  solutions  of  EDB  in  methanol  (several  microliters  into 
18  L  of  w'ater).  The  influent  water  was  stored  in  18  L  Teflon  bags 
to  prevent  loss  due  to  volatilization  and  covered  with  black  plastic 
to  prevent  photolysis  of  EDB.  Reduction  rates  of  EDB  in  spiked 
site  water  and  spiked  DI  water  were  identical. 

The  column  required  several  pore  volumes  before  the  EDB 
profile  reached  steady  state  conditions.  After  steady  state  was 
reached  (a  minimum  of  10  pore  volumes)  each  sample  port  was 
sampled  at  least  three  times  over  several  hours  (tens  of  pore  vol¬ 
umes).  The  flow  rate  was  then  increased  and  the  system  was 
allowed  to  reach  steady  state  again  before  sampling  was  resumed. 

Pilot  Scale  Packed-Bed  Reactor 

Fig.  2  shows  a  schematic  of  the  pilot  scale  packed-bed  reactor. 
The  influent  water,  FS-12  groundwater,  was  taken  from  the  9  in. 
(22.9  cm)  influent  pipe  at  the  existing  treatment  plant.  I'he  w^ater 
was  filtered  through  a  stainless  steel  strainer  (20x  20  mesh).  Two 
sets  of  sluiit>rf  valves  controlled  the  influent  and  effluent  flow,  one 
manual  the  other  a  normally  closed  solenoid  valve  controlled  by 
the  process  control  system.  A  rotometcr  indicated  flow  into  the 
system.  Flow  rale  was  controlled  by  a  NEMO  variable  speed, 


progressing  cavity  pump  (Netzsch  Inc.,  Exton,  Penn).  The  pump 
speed  was  controlled  with  an  AC  Tech  variable  speed  ac  motor 
drive  (AC  Tech).  The  lower  limit  of  the  NEMO  pump  was  1,000 
mL/min.  In  order  to  operate  the  system  at  lower  flow  rates  a 
portion  of  the  flow  was  shunted  through  a  bypass.  A  manual  valve 
controlled  the  volume  of  water  bypassed. 

The  emergency  solenoid  cutoff  valves,  the  pump  and  the  in¬ 
strumentation  were  controlled  using  a  process  control  software 
package,  Lookout  (National  Instalments,  Austin,  Tex.),  an  object- 
based  automation  software  system.  The  software  was  resident  on 
a  Pentium  PC  with  an  I/O  board,  AT-M10-16E-10  (National  In¬ 
struments,  Austin,  lex.)  and  digital  acquisition  software,  DAQ 
(National  Instruments,  Austin  Tex.)  A  pressure  transducer  (maxi¬ 
mum  pressure  1,000  psi,  accuracy ±0.25%,  Barksdale,  Los  An¬ 
geles),  and  a  stainless  steel  pressure  gauge  (maximum  pressure 
300±5psi,  McMaster-Carr,  Atlanta)  were  used  to  monitor  and 
display  influent  pressure.  The  flow  measured  by  the  flow  sensor 
(Flo-Sensor  102-7P,  100-2000  mL/min,  ±3%  FS,  McMillan  Co., 
Georgetown  Tex.)  was  read  by  the  process  control  system  and 
used  to  determine  the  pump  speed;  thus  the  flow  rate  remained 
constant  against  increasing  backpressure.  The  flow  was  manually 
checked  periodically.  The  pressure  drop  across  the  packed-bed 
reactor  was  measured  by  a  differential  pressure  sensor-transmitter 
(Orange  Research,  Milford  Conn).  If  the  differential  pressure  ex¬ 
ceeded  150  psi  the  pilot  plant  automatically  shut  off  and  an  emer¬ 
gency  pressure  relief  valve  opened.  Lookout  automatically  logged 
the  influent  pressure,  differential  pressure,  flow  through  the  flow 
transmitter,  and  the  calculated  flow  through  the  bypass  line. 

The  packed-bed  reactor  consisted  of  12  stainless  steel  canisters 
(length=53.3  cm.  diameter=8.9  cm,  empty  volumc  =  3,316  cnr\ 
Osmonics,  Minnetonka  Minn.)  plumbed  in  series  with  0.95  cm 
stainless  steel  tubing  and  Swaglok  fittings.  All  canisters  were 
identical  with  the  exception  of  the  first.  In  the  bottom  of  the 
canisters  were  stainless  steel  wire  mesh  screens  (32X32  mesh) 
followed  by  a  2. .5  cm  deep  layer  of  fine  sand  to  distribute  flow. 
The  rest  of  Itte  canister  was  filled  with  approximately  9.0  kg  of 
cast  iron  filings.  'The  exit  of  the  canister  was  packed  with  stainless 


JOURNAL  OF  ENVIRONMENTAL  ENGINEERING  /  JANUARY  2002  /  87 


Steel  wool  to  prevent  washout.  Samples  were  taken  after  each 
canister  and  at  the  influent  prior  to  the  flow  sensor.  A  check  valve 
after  the  last  canister  prevented  backflow.  The  system  effluent 
re-entered  the  9  in.  FS-12  treatment  system  influent  pipe  61  cm 
downstream  of  the  intake  point. 

The  first  canister  was  designed  to  remove  DO  with  minimal 
increase  in  pressure  drop.  It  had  a  2.5cmX38cm  stainless  steel 
wire  mesh  (32X  32  mesh  size)  center  tube  filled  with  fine  sand, 
surrounded  with  cast  iron  to  the  top  of  the  canister.  The  porosity 
of  the  sand  was  less  than  that  of  the  iron  filings,  so  preferential 
flow  was  through  the  iron.  As  the  iron  at  the  bottom  oxidized  and 
the  permeability  decreased,  the  water  flowed  up  the  center  post 
and  spread  into  the  fresh  iron.  Thus  as  a  portion  of  the  iron  was 
used  up  a  new  area  was  exposed.  Two  DO  scrubbers  were  in 
parallel  so  a  fresh  canister  could  be  put  on  line  without  interrupt¬ 
ing  the  flow.  Each  canister  was  used  until  the  pressure  drop 
reached  100  psi,  after  that  it  was  replaced.  The  useful  lifetime 
depended  on  flow  rate  and  influent  DO.  Canisters  2-12  were  kept 
in  service  throughout  the  duration  of  the  study.  The  DO  in  these 
canisters  never  exceeded  0.5  mg/L.  Sampling  was  not  performed 
for  at  least  24  h  after  switching  to  a  new  scrubber  canister  to 
allow  the  system  to  reach  steady  state. 

Pilot  plant  samples  were  collected  in  40  mL  amber  VOA  vials. 
Canisters  were  sampled  in  reverse  order,  12  to  influent.  Duplicate 
samples  were  taken  at  canister  6  and  the  influent.  The  vials  were 
completely  filled,  without  headspace,  with  no  preservative. 
Samples  were  shipped  overnight  with  icepacks  and  refrigerated 
until  extraction.  All  samples  were  extracted  as  described  below 
witliin  24  h  of  arriving  at  the  laboratory. 

Analytical 

Headspace  samples  were  collected  and  directly  injected  into  a 
Hewlett-Packard  5890  GC  equipped  with  a  FID  detector.  The 
column  was  a  DB-5  (30  m  i.d.  0.53  mm,  flow  rate  3.5  mL/min); 
the  oven  program  was  isothermal  at  40°C.  Relative  response  fac¬ 
tors  were  determined  for  EDB  (dimensionless  Henry's  constant  at 
20^C  was  0.035)  and  ethylene  (H=8.9)  (Sanders  1999),  by  add¬ 
ing  known  masses  to  160  ml  vials  containing  the  same  liquid/ 
vapor  ratio  as  the  reaction  vials.  Standards  were  equilibrated  at 
the  same  temperature  as  the  reaction  vials  (Rajagopal  and  Burris 
1999). 

The  method  of  sample  extraction  depended  on  the  expected 
concentration.  For  EDB  concentrations  less  than  20  ppb  (pilot 
scale  and  low  [EDB]  lab  scale  samples),  a  5  mL  aliquot  was 
transferred  to  an  8  mL  test  tube  and  '^0.50  g  of  NaCl  was  added. 
The  tube  was  capped  with  a  Teflon  lined  cap  and  mixed  to  dis¬ 
solve  the  salt.  The  aliquot  was  extracted  with  0.20  mL  of  THM 
grade  pentane  spiked  with  5.0  ixg/mL  DCB  as  an  internal  stan¬ 
dard,  and  shaken.  The  pentane  layer  was  transferred  to  a  100  p.L 
insert  in  a  2  mL  crimp  top  vial  for  analysis.  For  samples  of 
EDB>20ppb  (high  [EDB]  lab  scale  samples),  a  50  piL  sample 
was  immediately  spiked  into  autosampler  vials  containing  100  fxL 
10  ppm  DCB  as  an  internal  standard  in  acetonitrile.  At  higher 
concentrations  (10,000  ppb).  a  50  fxL  aliquot  was  immediately 
diluted  with  1  mL  of  acetoniirile  containing  10  ppm  DCB  in  a 
crimped  GC  vial. 

Samples  from  the  laboratory  columns  were  analyzed  for  EDB 
on  a  Hewlett-Packard  GC  equipped  with  an  electron  capture  de¬ 
tector  (ECD).  Samples  were  injected  in  splitless  mode  (split  after 
60  s)  at  150°C  onto  a  HP-1  column  (3  jxm  thickness,  0.53  mm 
i.d.,  15  m  length).  The  oven  temperature  program  was  40°C  for  2 
min  and  ramped  at  20^X7min  to  220X.  ECD  temperature  was 


Time  (hours) 

Fig.  3.  Reduction  of  ethylene  dibromide  (EDB)  by  ZVT  in  a  well- 
mixed  system.  [EDB]q=  1 1.8  p-M,  5  g  iron,  100  mL  DI,  headspace 
60  mL,  Vials  were  mixed  at  24  rpm  in  the  daik  at  20®C.  The  only 
organic  product  found  was  ethylene.  Control  (no  Fe^)  (•),  EDB  (O), 
C2H4  (A),  C  mass  balance  (A). _ 


280°C,  Helium  was  the  carrier  gas  (4.77  cm^/min)  and  the  detec¬ 
tor  make  up  gas  was  5%  methane  in  argon  (57  cmVmin).  The 
EDB  detection  limit  was  0.05  ppb. 

Dissolved  oxygen,  pH,  and  temperature  were  measured  in  the 
field  at  the  time  of  sample  collection.  Temperature  and  DO  were 
measured  using  a  YSI  model  58  Dissolved  Oxygen  Meter  with  a 
YSI  5905  BOD  probe  (YSI  Inc.  Yellow  Springs,  Ohio).  The  DO 
meter  was  calibrated  before  each  use  to  100%  air  saturation.  An 
Orion  250  A  electrode  was  used  for  pH  monitoring  (Orion  Re¬ 
search  Inc.,  Boston). 


Results  and  Discussion 
Batch  System 

Results  of  a  typical  experiment  in  the  batch  system  at  20°C  are 
given  in  Fig.  3.  The  only  organic  product  detected  was  ethylene 
and  the  carbon  mass  balance  was  near  100%.  Tliis  is  in  agreement 
with  Rajagopal  and  Burris  (1999)  who  studied  the  reduction  ot 
EDB  by  Fisher  40-mesh  cast  iron.  They  determined  that  ethylene 
was  the  major  product  of  EDB  reduction  and  that  no  bromoethane 
or  vinyl  bromide  was  formed. 

Pseudo-first-order  rate  constants  relative  to  surface  area 
were  determined  at  five  temperatures  (5,  10,  20,  30,  and  40®C). 
at  20°C  was  found  to  be  independent  of  initial  EDB  concentration 
(0.011-0.12  mM)  and  independent  of  mixing  rate  at  rotation 
speeds  greater  than  8  rpm  (data  not  shown).  The  portion  of  EDB 
in  the  headspace  was  less  than  l()7r  (Rajagopal  and  Burris  1999), 
thus  the  vapor-liquid  partitioning  was  not  expected  to  signifi¬ 
cantly  affect  the  reaction  kinetics.  An  Arrhenius  plot  [Eq.  (10)] 
returned  an  Arrhenius  parameter  /t  =  2.28Xl0^s  ^m  ^  and  an 
activation  energy  of  1^3—  55  ±  1 3  kJ/mol  (Rajagopal  and  Burris 
found  an  £:„=51  ±7  kJ/mol  for  l-ishcr  40-mesh  cast  iron). 

The  sorption  isotherm  for  EDB  was  measured  using  relatively 
high  concentrations  of  EDB  and  Miiall  amounts  of  metal  so  that 
the  system  had  time  to  equilibrate  and  samples  could  be  taken 
before  a  significant  fraction  of  EDB  was  reacted  (Fig.  4).  The 
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Fig.  4.  Sorption  isotherm  of  ethylene  dibromide  (EDB)  on  MB  cast 
iron  at  20''C.  The  line  was  fit  using  a  Langmuir  equation, 

=  29.53  and  EDB, ^,,=2. 17X1 0“^’  mol/nr. 


high  EDB  concentration  also  reduced  possible  interference  by 
reaction  products.  The  sorption  isotherm  was  fit  using  the  Lang¬ 
muir  equation  [Eq.  (4)]  EDB„^.,.^=2. 17X  10“^  mol/m“  and 
=  29.53  m^/mol.  Burris  et  ah  (1998)  showed  that  sorption  of  ha- 
logenated  organics  to  cast  iron  surfaces  occurs  at  both  reactive 
and  nonreactive  sorption  sites.  The  relative  concentration  of  reac¬ 
tive  sites  to  nonreactive  sites  on  the  cast  iron  was  not  determined. 
The  sorption  isotherm  in  Fig.  4  was  for  sorption  to  all  sites.  In  a 
well-mixed  system  mass  transfer  is  rapid  and  steady  state  condi¬ 
tions  arc  established  quickly.  It  was  assumed  that  the  change  in 
aqueous  EDB  concentration  was  due  to  reaction  at  active  sites 
and  sorption  from  the  aqueous  phase  to  the  open  active  site. 

Laboratory  Scale  Packed-Bed  Reactor 

Groundwater  received  from  FvS-12  MMR  was  filtered  and  spiked 
with  EDB  and  used  as  influent  solutions  for  the  lab  scale  studies. 
Pseudo-first-order  observed  removal  rate  constants  (/Cobs)  deter¬ 
mined  at  different  initial  EDB  concentrations  and  flow  rates  arc 
shown  in  Fig.  5.  The  rate  constants  were  calculated  after  steady- 
state  conditions  were  reached.  At  low  initial  EDB  concentrations 
values  increased  when  How  rate  increased.  At  higher  initial 
EDB  concentrations  values  were  lower  and  less  dependent  on 
flow  rates.  To  understand  the  effects  of  concentration  and  flow 
rate  on  observed  removal  rates  it  is  instructive  to  look  at  a  con¬ 
ceptual  model  of  heterogeneous  reactions. 

Model 

A  heterogeneous  reaction  can  he  thought  of  as  a  three  step  pro¬ 
cess,  (1)  adsoiption-dc.soi  piion  of  the  substrate  on  the  surface,  (2) 
reaction  of  the  surface  bound  substrate,  and  (3)  the  desorption  ol 
the  products  (Carberry  1976)  Assuming  steps  1  or  2  arc  rale 
limiting,  the  removal  rate  constant  of  an  aqueous  solute 
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Fig.  5.  Observed  removal  rate  constants  (A'obs)  of  EDB  in  lab  scale 
reactor  at  different  flow  rates  I,  5,  and  lO  mL/min,  and  initial  con¬ 
centrations,  6  (•),  1 00  (O),  1, 000  (▼),  10,000  ppb  (V). 


depends  on  mass  transfer  to  and  from  the  surface  and  the  reaction 
rate  on  the  surface  (intraparticle  diffusion  was  neglected). 

^  oh.s  ^ [  ^at]]  ^  [  E'surf]  (^^ 

where  (mVm“s)  =  rate  constant  of  mass  transfer,  ki 

(mVm^  s)  =  rate  constant  of  the  surtace  reaction,  a  =  ratio  of  solid 
surface  area  to  liquid  volume  (m^/m^),  and  C^m-ffmol/m^) 
=  concentration  of  solute  sorbed  to  the  solid  surface.  When  the 
rates  of  transfer  to  and  from  the  sorption  sites  are  equal  the  sys¬ 
tem  is  in  equilibrium  and  the  overall  mass  transfer  is  zero.  Arnold 
et  al.  (1999)  showed  that  mass  transfer  effects  were  greater  for 
substrates  with  faster  surface  reaction  rates.  Since  EDB  reacted 
relatively  rapidly  mass  transfer  was  expected  to  be  more  impor¬ 
tant  than  with  slower  reacting  species.  The  for  EDB  at  higher 
flow  rates  (5,  10  mL/min  in  Fig.  5)  indicated  that  mass  transfer 
did  affect  the  kinetics  and  that  this  effect  declined  as  aqueous 
concentration  increased.  In  order  to  elucidate  the  mass  transfer 
effects  in  packed-bed  reactors,  the  aqueous  EDB  concentration 
profiles  through  the  lab  scale  column  at  experimental  flow  rates 
and  concentrations  were  simulated  using  a  one  dimensional  non¬ 
equilibrium  model. 

Fig.  6  shows  the  aqueous  EDB  profiles  through  the  lab  scale 
reactor  at  two  flow  rates,  1  mL/min  and  10  mL/min  ([EDB]q 
=  6  ppb),  and  the  lines  that  were  fit  to  the  data  by  the  model 
described  below.  The  model  used  values  of  C,  S,  and  of 
1.15X10'^,  0.69,  and  12.7  (nrVnr),  respectively.  These  values 
were  determined  by  least  squares  regression  using  Powell’s  algo¬ 
rithm  (Micromath).  Average  values  of  C,  5,  and  were  calcu¬ 
lated  and  used  to  simulate  the  data  (solid  lines  in  Fig.  6).  The  r 
values  of  the  model  simulaticuis  of  (he  experimental  lab  scale  data 
were  >0.93. 

In  systems  with  high  liquit!  flow  rates  and  surface  reactions, 
partitioning  equilibrium  between  lic|uid  and  solid  phases  may  be 
retarded.  The  dependence  of  FDB  removal  in  the  lab  scale  col- 
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Fig,  6.  Modeling  of  lab  scale  packed-bed  reactor  data  at  two  flow 
rates  1  mL/inin  and  10  mL/min,  using  one-dimensional  nonequilib¬ 
rium  model.  The  model  parameters  C,  5,  and  were  determined 
by  fitting  the  model  to  the  data  sets  from  the  lab  scale  reactor. 
Duplicate  data  (O)  are  fit  by  the  model  (— ).  [EDB]o=6  ppb, 
C=1.15X10"^  5  =  0.69,  A:in=  12.7  mVml  _ 


umn  on  aqueous  solute  concentration  and  flow  rate  indicated  non- 
equilibrium  conditions.  Describing  these  nonequilibrium  systems 
required  a  kinetic  sorption  model.  A  bilinear  adsorption  model  is 
the  kinetic  version  of  the  Langmuir  isotherm  (Fetter  1993). 

=  -  ^forw[EDB](S^-  [EDB,„rf])  +  A:rev[EDB,^] 

(3) 

Mass  transfer  to  and  from  the  solid  surface  controls  the  change  in 
the  aqueous  concentration  of  EDB.  /cf^nv  is  the  rate  of  transfer  to 
the  surface  and  krev  is  the  rate  of  transfer  from  the  surface  to  the 
liquid,  EDBjjurf  is  die  concentration  of  EDB  on  the  surface  and 
is  the  EDBsurf  concentration  at  equilibrium.  The  can  be  found 
from  the  Langmuir  sorption  isotherm. 

—  ^ads[^i^^]  ... 

l  +  A:,ds[EDB] 

where  EDBn,ax==i'‘^^^ii^*J^  concentration  of  surface  sites,  (which 
is  approached  asymptotically  as  the  aqueous  concentration  in¬ 
creases)  and  A"ad.s=  equilibrium  constant  found  from  the  sorption 
isotherm  (Fig.  4). 

As  [EDBsurf]  approaches  the  equilibrium  surface  concentration 
the  first  term  on  the  right-hand  side  of  Eq.  (3)  approaches  zero 
and  the  rate  of  mass  transfer  declines.  At  higher  aqueous  phase 
concentrations  the  sorption  sites  on  the  surface  were  saturated. 
The  effects  of  higher  linear  velocities  on  removal  rates  became 
weaker  as  the  maximum  EDB  surface  concentration  was  ap¬ 
proached.  Thus  depending  on  the  aqueous  phase  concentration  the 
reaction  was  mass  transfer  limited,  surface  reaction  limited  or  a 
combination  of  the  two. 

The  second  driver  of  mass  transfer  was  the  flow  rate.  The  flow 
rate  (mVs)  was  transfonned  to  the  linear  velocity  (V),  which  was 
the  flow  divided  by  cross-sectional  area  of  the  column  (/4),  and 
porosity  (e). 


Table  1.  Linear  Velocities,  Reynolds  Numbers,  and  Reactor  Length/ 
Area  Ratios  (LI A)  for  Flow  Rates  in  Bench  Scale  Column  (1,5,  and 
10  mL/min),  Pilot  Plant  (300,  600,  and  2,000  mL/min),  and  Hypo¬ 
thetical  Full  Scale  Treatment  System  (2.3X  10^’  mL/min) _ 


Flow  (mL/min) 

V  (m/s) 

UA  (m"') 

1 

6.2X  10“^ 

0.025 

3.9 

5 

3.1X10"^ 

0.127 

3.9 

10 

6.2X  lO-* 

0.250 

3.9 

300 

1.5X  10“’ 

0.606 

10.3 

600 

3.0X  10"-’ 

1.21 

10.3 

2,000 

l.OX  10"- 

4.07 

10.3 

2,300,000 

7.3  X  10“^ 

3.0 

8.3 

Note:  The  model  can  be  scaled  up  as  long  as  the  L/A  and  N rc  are  within 


the  same  order  of  magnitude. 


The  particle  Reynolds  number  (A/^g)  was  dependent  on  the  linear 
velocity  (K),  the  average  cross  sectional  diameter  of  the  particles 
and  the  kinematic  viscosity  of  water  (nn) 


The  particle  Reynolds  numbers  for  flow  rates  used  in  this  study 
are  given  in  Table  1.  These  were  very  low  Reynolds  numbers  and 
indicated  a  laminar  flow  regime. 

Wilson  and  Geankopolis  (1964)  developed  an  empirical  equa¬ 
tion  for  mass  transfer  in  packed-bed  reactors.  Starting  from  this 
equation  the  value  of  was  estimated 

kt^r.=  CXNi,XSc^'^X^  (7) 


where  C  and  S  are  constants  whose  values  were  determined  by 
performing  a  least  squares  regression  of  the  column  data  for  the 
two  lowest  EDB  concentrations,  6  and  100  ppb,  at  three  flow 
rates.  Data  from  the  laboratory  scale  reactor  were  modeled  and 
the  values  of  C  and  S  were  varied  until  the  residual  differences 
between  the  experimental  EDB  data  and  the  predicted  EDB  con¬ 
centrations  reached  a  minimum.  The  6  and  100  ppb  data  showed 
the  greatest  dependence  on  mass  transfer  so  this  data  was  used  to 
fit  C  and  5,  obtaining  1.1 5X  10““^  and  0.69,  respectively. 

Dedb  is  the  diffusivity  of  EDB  in  water  at  a  specific  tempera¬ 
ture,  and  Sc  is  the  Schmidt  number.  The  Schmidt  number  is  the 
diffusivity  of  a  molecule  in  a  given  fluid. 


^forvv  is  dependent  on  system  parameters  such  as  particle  diam¬ 
eter  and  porosity,  which  were  constant  and  variables  dependent 
on  temperature  (Sc,  t),  Dedb)  and  linear  velocity  (A/Rg),  In  the  lab 
scale  studies  the  reaction  temperatures  (22“23°C)  are  higher  than 
in  the  pilot  scale  system  (9-12‘'C).  The  temperature  dependence 
is  included  so  the  model  could  be  applied  to  both  sets  of  data. 

When  the  surface  concentrations  were  low  /:forw[EDB](5sat 
-[EDBsijrt])>^rev[EE>Bsurf]-  Under  these  nonequilibrium  condi¬ 
tions  the  sorption  of  EDB  to  the  surface  was  greater  than  desorp¬ 
tion.  This  created  transient  conditions  where  mass  transfer  limi¬ 
tations  were  diminished  and  removal  rates  much  higher  than  those 
observed  at  higher  aqueous  phase  concentrations  or  lower  linear 
velocities  were  achieved.  Such  transient  conditions  were  de¬ 
scribed  by  Fetter  (1993)  for  the  sorption  of  inorganic  solutes  onto 
mineral  surfaces.  On  reactive  surfaces  where  surface  solute  con- 
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centralions  were  affected  by  a  chemical  process  (k^)  and  a  physi¬ 
cal  process  these  transient  conditions  are  maintained  longer. 

At  equilibrium,  In  the  model  it  was  assumed 

that  ^,ev=^forw/^ads-  Bccausc  the  lab  scale  data  inputted  into  the 
model  was  collected  after  the  system  reached  steady  state,  it  was 
necessary  to  estimate  the  concentration  of  EDB  on  the  surface 
in  the  initial  model  parameters.  The  boundary  conditions  were 
0<[EDB,u,.f]o<5j.j^t.  ^rev  was  not  a  true  constant  as  written,  to 
account  for  this  [EDB^^Jo  was  set  to 


^  EDB,,,/^i,rEDB] 
[EDB,,,fJo-  ,+^jedB] 


(9) 


where  A'in=  initial  (nonequilibrium)  ratio  of  /cfonv/^rev  This  value 
was  found  by  applying  least-squares  regression  to  the  lab  scale 
packed-bed  reactor  data.  A  value  of  12.7  minimized  the 
error  between  the  observed  data  and  the  concentrations  predicted 
by  the  model.  After  the  first  iteration,  the  value  of  EDBs^„f  was 
calculated  using  Eq.  (11). 

The  concentration  of  EDB  on  the  reactive  surface  was  depen¬ 
dent  on  mass  transfer  and  the  surface  reaction.  The  rate  of  surface 
reaction  was  calculated  using  and  A  from  the  batch  experi¬ 
ments 


fc,=/1exp|^|  (10) 

Combining  the  mass  transfer  terms  and  the  surface  reaction 
term  gave  the  rate  of  change  of  [EDB^^^rf]- 

r/fEDB^fl 

j/  =  ^fonv.^[EDB,q]  -  ^,ev[EDB,,pf]  -  it  i^[EDB,,,f] 

(11) 

^^ame  from  Eq.  (7).  Eq.  (10)  gave  the  value  for  k^,  and 
^rcv^^forw^^jids-  Thc  tit  of  the  lab  scale  data  was  found  to  be 
independent  of  k]  indicating  that  mass  transfer  controlled  the  re¬ 
moval  rate  of  EDB  in  the  lab  scale  experiments. 

Pilot  Scale  Packed-Bed  Reactor 

After  a  brief  period  of  in-house  testing,  the  pilot  scale  system  was 
transported  to  MMR  in  October  1998.  The  pilot  plant  was  in  use 
for  8  months  from  October  1998  to  June  1999,  more  than  82,000 
L  of  groundwater  were  treated.  The  EDB  influent  concentration 
declined  throughout  the  length  of  the  experiment  due  to  ongoing 
treatment  of  the  plume,  [EDB]o”7-3  ppb.  Dissolved  oxygen 
was  removed  from  the  innuent  (DO<0.5mg/L)  in  the  first  can¬ 
ister.  The  pH  increased  from  initial  values  of  6. 8 ±0.5  to  a  steady 
value  of  7.3 ±0.5  after  the  first  canister  and  remained  relatively 
constant  throughout  the  system  at  all  flow'  rates. 

Three  influent  How  rates  were  used;  300,  600,  and  2,000  mL/ 
min.  Many  pore  volumes  (100-400)  were  required  for  the  con¬ 
centrations  to  reach  steady  state.  TJie  number  of  pore  volumes 
required  depended  on  flow’  rate  and  the  concentration  of  EDB. 
Sorption  to  nonreactivc  surfaces  has  been  showm  to  cause  tlie 
retardation  of  solutes  in  porous  media,  however,  once  steady  slate 
conditions  are  reached,  nonreactive  sorption  will  not  affect  the 
extent  of  reaction  (Hatfield  cl  al.  1996). 

Fig.  7  show's  EOB  concentrations  in  the  pilot  system  for  a  flow 
of  600  mL/min.  The  data  from  37  to  645  pore  volumes  show's  the 
delay  in  reaching  steady  slate  conditions.  Once  steady  state  wais 
reached  the  EDB  prolile  through  the  pilot  plant  was  stable  ha 
hundreds  of  pore  volumes.  The  target  effluent  concentration  ol 
[EDB]  =  0.02  ppb  w'as  not  reached,  'fhe  removal  efficiencies  for 


Fig.  7.  Ethylene  di bromide  (EDB)  concentration  profile  through  the 
pilot  plant,  y=  600mL/inin.  The  concentrations  at  each  sample  point 
increase  over  time  (from  37  to  1943  pore  volumes)  until  steady  state 
conditions  are  reached.  Once  steady  state  is  reached  the  profile  is 
stable  for  hundreds  of  pore  volumes,  645-1943  pore  volumes.  The 
initial  rapid  drop  is  due  to  nonequilibrium  mass  transfer.  The  removal 
of  EDB  from  canisters  2-12  is  dominated  by  the  rate  of  the  surface 
reaction. 


flows  of  300  mL/min,  600,  and  2,000  after  the  last  canister  (12), 
were  73,  49,  and  39%  respectively. 

The  concentration  profile  through  the  pilot  plant  showed  tw'O 
regions.  The  first  was  a  sharp  reduction  in  EDB  concentration 
after  the  first  canister  (at  300  mL/min,  after  the  second  canister). 
After  that,  the  EDB  concentrations  declined  slowly  and  steadily. 
Pseudo-first-order  rate  constants  for  EDB  removal  in  these  two 
regions  show'  that  in  the  first  the  rate  constants  increased 

w'ith  increasing  flow  rate.  In  the  second  region  the  rate  constants 
(^sec)  independent  of  flow'  rate  (Table  2  and  Fig.  7). 

This  dependence  on  linear  velocity  in  the  first  region  was  simi¬ 
lar  to  that  observed  in  the  lab  scale  studies.  Tliis  showed  that  mass 
transfer  dominated  in  the  early  stages  of  the  pilot  scale  reactor  but 
once  the  transient  nonequilibrium  conditions  had  passed  the  rate 
of  surface  reaction  controlled  the  kinetics. 

The  data  from  the  pilot  scale  experiments  were  simulated 
using  the  model  developed  from  the  batch  isotherm  and  the  mass 
transfer  equations  from  the  lab  scale  studies.  The  pilot  scale  sys¬ 
tem  had  a  different  reactor  configuration  (Table  1),  lower  process 
temperatures,  and  lower  initial  EDB  concentrations  than  those 


Table  2.  Pscudo-First-Order  Rate  Con.stants  From  Pilot  Plant  Experi¬ 
ments 


Flow  rate  (L/min) 

V  (cm/min) 

Wr. 

^■|ir.u  (inin  ') 

(min  ') 

0.3 

4.8 

0.42 

0.052±.004 

0.014±.004 

0.6 

9.6 

0,84 

0. 1 38 ±.048 

0.0 13  ±.005 

2.0 

32.2 

2.79 

().249±.086 

0.019±.012 

N'oie:  A'lirsi  rate  constants  were  calcufucd  using  the  first  three  data  points 
(induent  to  canister  2).  constants  were  derived  from  the  steady  state 
uihylcne  dibromide  concentrations  in  canisters  2-12. 
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Fig.  8.  Pilot  scale  packed-bed  reactor  data  and  model  simulations. 
K=300  mL/min  (data  M,  model  — ^),  [EDB]o  =  22.3  fxM,  T 
^ll.O'^C,  r^  =  0.99;  K=600  mL/min  (data  A,  model  — ),  [EDB]o 
=  15.4  fxM,  r=]1.3^C,  r^  =  0.98;  V=  2,000  mL/min  (data  ♦, 
model  — [EDB]o=  15.4  fxM,  r=12.2°C,  r^  =  0.93. _ 


Lab  scale  column  studies  varying  flow  rates  and  initial  EDB 
concentrations  exhibited  characteristics  of  nonequilibrium  sorp¬ 
tion.  Observed  removal  rates  of  EDB  in  aqueous  solutions  in¬ 
creased  at  higher  flow  rates.  Increasing  [EDB]o  decreased  the 
apparent  removal  rates  (Fig.  5).  Modeling  of  the  system  indicated 
that  under  nonequilibrium  conditions  there  was  rapid  sorption  to 
surface  sites,  at  higher  [EDBJo  and  longer  retention  times  (lower 
V),  this  phenomenon  diminished. 

The  pilot  plant  experiments  were  run  for  several  months  at 
three  flow  rates.  The  longer  hydraulic  retention  times  in  the  pilot 
plant  relative  to  the  bench  scale  columns  revealed  that  alter  the 
transitory  nonequilibrium  conditions  removal  was  dominated  by 
the  rate  of  the  surface  reaction.  These  transitory  conditions  com¬ 
plicated  the  analysis  of  the  system  kinetics.  Studies  done  with 
short  columns  over  predicted  removal  rates  and  would  have  lead 
to  under  sizing  of  packed-bed  reactors. 

Packed-bed  ZVI  reactors  may  provide  an  alternative  to  the 
current  best  available  technology,  granular  activated  carbon 
(GAC).  The  advantages  of  ZVI  included  the  detoxification  of 
EDB  and  lower  operation  and  maintenance  costs  the  GAC. 
Packed-bed  ZVI  reactors  may  be  useful  treatment  options  for 
pump-and-treat  remediation  of  EDB  contaminated  groundwaters. 
Jetting  or  hydraulic  fracture  placement  of  ZVI  PRB  may  also 
provide  in  situ  treatment  options  for  EDB  contaminated  ground- 
water. 


used  in  the  lab  scale  experiments.  Variables  such  as  [EDB]o, 
temperature,  and  flow  rate  were  entered  into  the  model  as  initial 
parameters.  The  value  was  fitted  by  the  least  squares  method 
using  data  from  all  three  flow  rates  and  was  determined  to  be  22 
(m^/m^).  Fig.  8  shows  the  experimental  data  and  the  model  pre¬ 
dictions  of  three  data  sets  (flow=300,  600,  and  2,000  mL/min). 
The  correlation  between  data  and  model  were  strong  (r^ 
>0.93).  The  model  was  able  to  predict  the  two  regions  in  the 
EDB  profile  that  were  seen  in  the  pilot  plant  data. 

The  model  was  used  to  design  a  hypothetical  above  ground 
treatment  system  for  conditions  similar  to  FS-12  MMR.  The  sys¬ 
tem  was  scaled  up  by  keeping  and  the  reactor  length/cross 
sectional  area  ratio  (17 A)  within  an  order  of  magnitude  of  those  of 
the  lab  and  pilot  scale  reactors  (Table  1).  Using  the  current  influ¬ 
ent  conditions  of  FS-12  MMR,  [EDB]o=3  ppb,  and  flow 
rate=612  gal/min  (2.3  mVs),  a  packed-bed  reactor  with  a  cross 
sectional  area  of  3.0  m,  a  length  of  1 10  m,  and  a  hydraulic  reten¬ 
tion  time  of  3. 1  hours  would  be  required  to  reach  an  effluent  EDB 
concentration  less  than  20  ppt. 

A  hypothetical  hydraulic  fracture  placement  PRB  was  also 
modeled.  The  average  groundwater  flow  rate  at  FS-12  MMR  is 
1.5  ft/day  (5.3  |jl/s)  (MMR  2000).  Hocking  et  al.  (1998)  reported 
that  the  average  thickness  of  a  PRB  emplaced  by  hydraulic  frac¬ 
ture  was  0.089  m  (3.5  in.).  Assuming  a  groundwater  temperature 
of  irC,  and  that  e  and  a  were  the  same  as  the  packed-bed  reac¬ 
tor,  an  initial  EDB  concentration  of  3  ppb,  a  wall  thickness  of 
0.164  m  (6.5  in.)  would  be  required  to  achieve  the  EDB  MCL 
level  of  20  ppt. 


Conclusions 

Mass  transfer  and  .sorption  played  significant  roles  in  the  kinetics 
of  EDB  removal  by  cast  iron  packed-bed  reactors.  Initially  there 
was  significant  nonreactive  sorption  of  EDB  to  the  iron  surface 
(Fig.  7).  Tens  of  pore  volumes  were  required  to  reacli  partitioning 
equilibrium. 

92 /JOURNAL  OF  ENVIRONMENTAL  ENGINEERING  /  JANUARY  2002 


Acknowledgments 

Funding  was  provided  by  the  Air  Force  Center  tor  Environmental 
Excellence  and  the  Air  Force  Office  of  Scientific  Research.  The 
writers  would  also  like  to  thank  Maria  Zerwas,  Paul  Moody,  and 
Robert  Nichols  of  Applied  Research  Associates,  Inc.  for  their  in¬ 
puts. 

Notation 

The  following  symbols  were  used  in  this  paper. 

A  =  cross-sectional  area  of  reactor  (m^); 

A  =  Arrhenius  parameter; 
a  =  iron  surface  area/liquid  volume  (m“/m'^); 

C  —  mass  transfer  constant; 

Cgurf  =  concentration  of  sorbed  solute  (mol/m"); 

Dedb  ^  diffusivity  of  EDB  in  water  (m/s); 

-  particle  average  cross  sectional  diameter  (m); 

Ea  =  activation  energy,  kJ/mol; 

EDBf^ax  “  maximum  concentration  of  surface  sites 
(moiym“);  ^ 

EDBsurf  —  concentration  of  sorbed  EDB  (mol/m“); 

/  =  flow  rate  (mVs); 

^ads  ”  equilibrium  constant  (mVm^); 

=  initial  ratio  of rev » 

^first  ”  observed  rate  constant  of  EDB  removal  in  first 
region  (min“'); 

^forw  “  sorption  rate  constant  (s“^); 

=  general  mass  transfer  rate  constant  (mVnrs); 

^obs  ~  observed  removal  rate  constant  of  lab  scale  cx- 
perimcnls  (s'  ^); 

-  desorption  rate  constant  (s“*); 

A^sec  ^  observed  rate  constant  of  EDB  removal  in  sec¬ 
ond  region  (min"’); 

ki  -  rale  consiani  of  surface  reaction  (m^/m“s); 

L  =  reactor  length  (m); 


=  particle  Reynolds  number; 

S  =  Reynolds  number  exponent; 

at  equilibrium  (mol/m~); 
Sc  =  Schmidt  number; 

V  —  linear  velocity  (m/s); 
e  -  porosity;  and 

T)  =  kinematic  viscosity  of  water  (m7s). 
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